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The charged-particle multiplicity distribution and its moments have been mea- 
sured, for all hadronic as well as for light-quark and b-quark events in e+e — col- 
lisions at the Z mass. The H q moments derived from the charged-particle multi- 
plicity distribution are known to exhibit quasi-oscillations when plotted versus the 
order of the moment. This behavior is predicted by the NNLLA of perturbative 
QCD for the parton level and, under the assumption of LPHD, also for the hadron 
level. Using the jet multiplicity distributions in order to vary the dependence on 
the LPHD hypothesis, we find, however, that at our energy the oscillations only 
appear for non-perturbative scales. In the absence of confirmation of pQCD, we in- 
vestigate a more phenomenological answer in the possibility that the features seen 
in the H q behavior could be due to the fact that the charged-particle multiplicity 
derives from a superposition of final states related to the topology of the events. 
Therefore, the analysis is repeated using charged-particle multiplicity distributions 
originating from 2-jet and 3-jet events for the full, light- and b-quark samples. 



1 Introduction 

Although the number of charged particles is only a global measure of the 
characteristics of the final state of a high-energy collision, it has proved a 
fundamental tool in the study of particle production. Independent emission 
of single particles leads to a Poissonian multiplicity distribution. Deviations 
from this shape, therefore, reveal correlations a. The shape of the charged- 
particle multiplicity distribution analysed with the ratio of cumulant factorial 
moments to factorial momentsB, H q , is known to reveal quasi-oscillations 13, 
when plotted versus the order q, with a first minimum at q = 5. 

The usual way to interpret this result is to refer to perturbative QCD, 
which provides us with calculations for the H q of the parton multiplicity 
distribution □. The Next to Next to Leading Logarithm Approximation 
(NNLLA) , which has the most accurate treatement of energy-momentum con- 
servation, predicts for the H q a negative first minimum near 5 followed by 
quasi-oscillations. This behavior may be expected for the charged-particle 
multiplicity distribution under the Local Parton-Hadron Duality (LPHD) hy- 
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pothesis, which assumes that the hadronization does not distort the shape of 
the multiplicity distribution. 

However, this result can also be interpreted in a more phenomenological 
way by viewing the shape of the charged-particle multiplicity distribution as 
a superposition of different types of event like 2-jet and 3-jet eventsO. This 
can be investigated using rather simple parametrizations, as a weighted sum 
of two Negative Binomial Distributions (NBD), each NBD carrying parame- 
ters (mean, n and dispersion, D) taken from the experimental 2-jet or 3-jet 
charged-particle multiplicity distributions and using as relative weight the 2- 
jet fraction (a 2 j t), 2NBD ful i samp i c = a2jctNBD 2jc t + (1 - a 2 jct)NBD3j ot . A 
similar parametrization can also be tested using light- and b-quark events!, 
instead of 2-jet and 3-jet events. 

The test of the two approaches is done by measuring charged-particle 
multiplicity distributions and their moments for the full, light- and b-quark 
samples. These samples were also subdivided into 2-jet and 3-jet events ob- 
tained from various y cu t values. _ 

This analysis is based on data collected by the L3 detectorLl in 1994 
and 1995 at the energy of the Z. The data sample corresponds to approx- 
imately one million selected hadronic events. A b-tag_algorithm is used to 
discriminate between light-(udsc) and b-quark eventsa. Furthermore, the 
resulting multiplicity distributions are fully corrected for selection, detector 
inefficienciesBll3li3 and light- or b-quark piLrity. It has been shown that the H q 
moments are very sensitive to truncationlii Since we want to compare a large 
variety of multiplicity distributions, we have to make sure that all distribu- 
tions are affected by the truncation in the same way. Therefore, the truncation 
is defined as the fraction of events removed in the tail of the charged-particle 
multiplicity distribution of the full sample (events with multiplicity larger 
than 48 are removed), and the same fraction of events is then removed in all 
charged-particle multiplicity distributions studied. 

2 Test of the pQCD approach 

The H q for the charged-particle multiplicity distribution of the full, light- and 
b-quark samples (figure exhibit a first negative minimum at q — 5 and quasi- 
oscillation for greater q. The H q measured from the light-quark sample are 
found to agree very well with those of the full sample, while slight differences 
exist, mainly at low q, for the H q measured from the b-quark sample. The 
observed behavior is similar to that predicted by the NNLLA. However, also 
JETSET EJ agrees very well with all the data samples, even though the parton 
shower of JETSET does not use NNLLA. The same behaviour is found when 
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Figure 1. H q of the charged particle mul- Figure 2. H q for various Monte Carlo op- 
tiplicity distribution. tions 

using other parton generation and fragmentation models and even when we use 
matrix element production of qq only, even with independent fragmentation. 
This shows us that the H q behavior can be reproduced without the need for 
the NNLLA of pQCD. 

Furthermore, our analysis of jet multiplicity obtained at perturbative en- 
ergy scales (> 1 GeV), where pQCD predictions for the behaviour of H q should 
be directly accessible, did not show any of the pQCD predictions made for 
the H q E£l . Therefore, the H q behavior seen for the charged-particle and 
jet multiplicity distribution at non-perturbative energy scales only, appears 
unrelated to the behavior of the H q calculated in NNLLA. 



3 The phenomenological approach 

This approach relies on the assumption that we can view the charged-particle 
multiplicity distribution as a superposition of distributions originating from 
various processes related to the topology of the event, as 2-jet, 3-jet, light- 
or heavy-quark events. Assuming that each of these processes can by itself 
be described by a relatively simple parametrization as the NBD, the charged- 
particle multiplicity distribution of the full sample would then be a weighted 
sum of all the contributions. All together, these various contributions would 
explain the shape of the charged-particle multiplicity distribution and, hence, 
the H q behavior. We checked mainly two hypotheses with this parametriza- 
tion. 
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Figure 3. Charged-particle multiplicity Figure 4. H q of the charged-particle mul- 

distributions for b-quark events and 2-jet tiplicity distribution for b-quark with 

and 3-jet b-quark events, together with the H q calculated from the 2NBD 

their parametrizations. parametrizations 

1. The first assumes that the shape of the charged-particle multiplicity 
distribution of the full sample arises from the superposition of 2-jet and 3- 
jet events— Our 2-jet and 3-jet samples were obtained using the Durham 
algorithmic for a set of six y cut values. As parameters for the NBD's, we 
used the means and dispersions calculated from the experimental 2-jet and 
3-jet charged-particle multiplicity distributions. The relative weight between 
the two NBD's was taken to be the fraction of 2-jet events for a given y C ut 
value. This gives us a fully constrained 2NBD parametrization of the full 
sample. The resulting \ 2 are given in the left half of table 1. Since the H q 
moments from charged-particle multiplicity distributions of full, light- and 
b-quark samples are very similar, we also tested this hypothesis on light and 
b-quark samples separately, isolating in these cases the 2-jet and 3-jet events 
from the light- and b-quark samples. The resulting x 2 f° r the b-quark sample 
are given in the right half of table 1. We find amazingly good \ 2 f° r the 
2NBD parametrizations of the full (column 2 of table 1), light- (not shown) 
and b-quark (column 5 of table 1) samples, (see also figure || for the case of 
the b-quark sample). We calculated the H q from the parametrizations and 
also these are found to be in good agreement with the H q measured for the 
full, light and b-quark samples (figure ^ for the b-quark sample). However, 
none of the NBD parametrizations are able to describe any of the individual 
2- or 3-jet charged-particle multiplicity distributions themselves, even though 
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Figure 5. H q of the charged-particle mul- Figure 6. Evolution of H-j as a function of 

tiplicity distribution, for all events and for y cut for 2-jet and 3-jet, light- and b-quark 

2-jet and 3-jet events (y cut = 0.015) events 

the x 2 is seen to decrease when the purity increases. 2. We also attempted 
to parametrize the 2-jet and 3-jet charged-particle multiplicity distributions, 
and as a consistency check the full sample, by a superposition of light and 
b-quark events, using in that case.^ relative weight between the two NBD's, 
the fraction of b-quark events, i?6tll. Results are summarized in table 2. We 
don't find any agreement at all, neither for 2NBD parametrization of the 
full sample which has a x 2 /dof near 13. This constitutes by its failure a 
good check of the method, since it shows that not all combinations of two 
NBD's agree with the data. We extended the study of the 2-jet and 3-jet 
samples to the measurement of their H q moments. We find that, even if the 
oscillations are still there, their amplitudes are far smaller than for the full 
sample (figure |^) . Furthermore the size of the oscillations decreases when the 
purity in 2 jet (3 jet) in the 2-jet (3-jet) sample increases. Differences at low q 
(mainly for q < 8) are found between the H q of 2-jet and 3-jet events, and also 
large differences are seen for fixed q when the Hj moment is plotted versus 
2/cut (figure^). We see differences between H q of 2-jet and 3-jet events of the 
light- and b-quark samples, however the oscillations are comparable to those 
of 2-jet or 3-jet events of the full samples. All together, this supports the 
phcnomenological approach when we assume that the shape of the charged- 
particle multiplicity distribution arises from a superposition of 2-jet and 3-jet 
events. 
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Table 1. \ 2 between the 2-jet, 3-jet parametrization and their experimental counterpart 





X^/dof for the full sample 


xVdof for the b-quark sample 


J/cut 


2NBD a n 


NBD 2 jot 


NBD 3jct 


2NBD all 


NBD 2jot 


NBD 3jct 


0.03 


1.4 


56 


8.5 


2.6 


19 


2.2 


0.015 


0.6 


36 


15 


1.4 


13 


3.4 


0.01 


20.4 
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6.3 


0.93 


9.8 


4.4 


0.006 


0.7 


16 


36 


0.53 


5.5 


6.4 


0.004 


1.5 


9.4 


57 


0.53 


3.5 


8.5 


0.002 


5. 


5.5 


119 


3 


0.8 


15 



Table 2. \ 2 between the light-, b-quark parametrization and their experimental counterpart 





xVdof 


J/cut 


2NBD 2 jot 


2NBD 3 j Ct 


0.03 


74 


8 


0.015 


46 


14 


0.01 


39 


9 


0.006 


23 


33 


0.004 


14 


48 


0.002 


5 


100 



4 Conclusions 

The oscillatory behavior of the H q moments of the charged-partiele multi- 
plicity distribution is usually interpreted as a confirmation of NNLLA, but 
investigations performed on different models of parton generation and for dif- 
ferent fragmentation models have shown similar oscillatory behavior in all 
cases. Furthermore, the analysis of the H q of the jet multiplicity distributions 
reveals that this behavior appears only for very small y cut , corresponding to 
energy scales < 100 MeV, far from the perturbative region. This gives us 
strong indications that the oscillatory behavior is not related to the behavior 
predicted by the NNLLA. In search of an alternative origin of this H q behav- 
ior we have, therefore, investigated a more phenomenological answer which 
assumes that the shape of the multiplicity distribution results from the su- 
perposition of 2-jet and 3-jet events. Using a weighted sum of 2 NBD's as a 
parametrization, we found very good agreement for both the charged-partiele 
multiplicity distribution and its H q moment. This supports the idea that the 
main feature in the shape of the charged-partiele multiplicity still visible in 
the final states are due to the presence of hard gluon radiation and to the 
hadronization. 
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